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We expand the most general lattice Dirac operator D in 
a basis of simple operators. The Ginsparg- Wilson equation 
turns into a system of coupled quadratic equations for the ex- 
pansion coefficients. Our expansion of D allows for a natural 
cutoff and the remaining quadratic equations can be solved 
numerically. The procedure allows to find Dirac operators 
which obey the Ginpsparg- Wilson equation with arbitrary 
precision. 
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In the last two years we have witnessed exciting devel- 
opments for chiral fermions on the lattice (see Q for some 
reviews). It was found that the perfect action approach 
H leads to the fixed point Dirac operator || which obeys 
the Ginsparg- Wilson equation Q. Also the overlap for- 
malism for chiral fermions on the lattice Q has led to an 
explicit expression for a Dirac operator which solves the 
Ginsparg- Wilson equation, the so-called overlap operator 
P). It was realized, that based on the Ginsparg- Wilson 
relation it is possible to redefine the chiral symmetry for 
lattice fermions and this new construction has spurred a 
wealth of interesting work Hj . 

Unfortunately the numerical implementation of the 
two known solutions of the Ginsparg- Wilson equation 
seems to be rather troublesome. Computing the inverse 
square root which appears in the overlap operator is a 
numerically challenging and expensive task. For the per- 
fect action the situation is even worse, since already the 
blocking procedure necessary for its computation is ex- 
tremely costly and so far has been implemented only in 
two dimensions 0. 

In this letter we present a new approach to the 
Ginsparg- Wilson equation: The first step is to expand 
the most general Dirac operator D in a basis of simple op- 
erators on the lattice. In the second step we insert the ex- 
panded D into the Ginsparg- Wilson equation which then 
turns into a system of coupled quadratic equations for the 
expansion coefficients of D. Solving the Ginsparg- Wilson 
equation is equivalent to finding solutions for this system 
of quadratic equations. Any solution of the Ginsparg- 
Wilson equation, such as the overlap operator or the 
perfect action, corresponds to a solution of our system 
of quadratic equations. Many interesting aspects can be 
addressed in this new framework, such as e.g. the ques- 
tion if there is a continuous manifold of solutions in the 
space of coefficients of D. 



Here we concentrate on exploring the possibilities for 
solving our system of quadratic equations numerically 
and using this approach for constructing new solutions 
of the Ginsparg- Wilson equation. We show that our ex- 
pansion of D provides for a natural cutoff which turns the 
quadratic equations into a simple finite system. Solving 
this system numerically gives the values for the expansion 
coefficients. This procedure allows to construct approxi- 
mate solutions which obey the Ginsparg- Wilson equation 
with arbitrary precision. 

Before we outline the details of our construction, let us 
begin with some remarks: Usually the derivative term on 
the lattice is discretized by the following nearest neighbor 
term (U^x) G SU(N) and we set the lattice spacing to 
1): 

1 4 

o [^W^+fcs ~ Un(x-fi)- l 6 x -n,y] . (1) 

However, it is perfectly compatible with all the symme- 
tries to instead discretize the derivative term using 
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U^{x-(iy 1 U^x~2fi)- 1 6 x - 2 fi,y , (2) 



and there are many more terms that are eligible. Thus an 
ansatz for the most general Dirac operator D must allow 
for a superposition of all of the possible discretizations 
for the derivative term. In order to remove the doublers 
we also have to allow for terms that come with 1 in spinor 
space and again we allow for all possible terms. We gen- 
eralize our D further, by including all terms also for the 
remaining elements T a of the Clifford algebra, i.e. we in- 
clude also tensor, pseudovector and pseudoscalar terms. 
From our examples ([j]) , @ it can be seen that the terms 
are characterized by the orientations of the link variables 
U )Jl (x). The ensemble of links supporting the gauge trans- 
porters can be viewed as paths and we can denote our 
ansatz for the most general lattice Dirac operator in the 
following form: 
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To each generator r Q of the Clifford algebra and to each 
pair of points x, y on the lattice we assign a set V x<y of 
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paths p connecting the two points. Each path is weighted 
with some complex weight cp* and the construction is 
made gauge invariant by including the ordered product 
of the gauge transporters Ui G SU(N) for all links I of p. 
The action is then given by S = J2 X y 1 Px^x,y' l Pyi where 
x and y run over all of the lattice. 

The next step is to impose on S the symmetries which 
we want to maintain: Translation and rotation invariancc 
and invariance under C and P. In addition we require 
our D to be 75-hermitian, i.e. we require 75-D75 = D* . 
This property can be seen to correspond to what leads 
to the CTP theorem in Minkowski space, i.e. the vector 
generators 7^ come with a derivative term etc. (compare 
Equation (0) below). 

Translation invariance requires the sets V" y of paths 
and their coefficients to be invariant under simultane- 
ous shifts of x and y, while rotation invariance implies 
that a term and its rotated image have the same weight. 
Parity implies that for each path p (with coefficient c^) 
we must include the parity-reflected copy with coeffi- 
cient Sp arity ■ Cp where the signs Sp arity are defined by 

74 la 74 = s parity ' "a- 

More interesting are the symmetries C and 75- 
hermiticity. It is easy to see that both of them imply 
a relation between the coefficient for a path p and the 
coefficient of the inverse path p . Implementing both 
these symmetries firstly restricts all coefficients to be 
either real or purely imaginary. Secondly we find that 
the coefficient for a path p and the coefficient for its 
inverse p~ x are equal up to a sign s" harge defined by 
CY a C = s" harge ■ rj, where T denotes transposition and 
C is the charge conjugation matrix. To be specific, we 
use the chiral representation for the 7^ and C = i^li- 
When implementing all these symmetries we find that 
paths in our ansatz become grouped together where, up 
to sign factors, all paths in a group come with the same 
coefficient. 

Before we present the explicit form of the most general 
Dirac operator with the above symmetries, let us first in- 
troduce an abbreviated notation: The terms in the Dirac 
operator are determined by the corresponding clement 
of the Clifford algebra and the paths with their relative 
signs. For the two examples (Eh, (0) the paths are given 
by one (two) steps in positive /i-direction and one (two) 
steps in negative //-direction with a relative minus sign. 
We now denote an arbitrary path of length n on the lat- 
tice by < ii, Z2, --In >j with the k giving the directions 
of the subsequent hops, i.e. k G {±1, ±2, ±3, ±4}. So 
for example the terms occurring in (0) are denoted by 
< /i, /i > and — < — (j,, —fj, >, where the minus in front of 
the second path is due to the relative minus sign in (0). 

Using this notation we can now write the most general 
Dirac operator on the lattice in the form: 
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The variables U run over all of {±1, ±2, ±3, ±4} unless 
restricted otherwise. We use the abbreviation s(/j) = 
sign(Zj). By e we denote the totally anti-symmetric ten- 
sors with 2,3 and 4 indices. We choose the normaliza- 
tion of the elements of the Clifford algebra such that the 
elements appear as all possible products of the 7^ with- 
out any extra factors of i. For this normalization the 
symmetries C and 75-hermiticity render all coefficients 
Si, Vi, ti, a,i and pi real. 

The above mentioned structure of paths appearing in 
groups is obvious in (|j) and we ordered the terms accord- 
ing to the length of the paths they contain. The paths 
in each group are related by symmetries and up to the 
sign factors have to come with the same real coefficient. 
All paths within a group have the same length. It has to 
be stressed, that in (Q) for each element of the Clifford 
algebra we show only the leading terms, i.e. the terms 
with the shortest paths, of an infinite series of groups of 
paths. The dots indicate that we omitted groups with 
paths that are longer than the terms we show. 

Let us now insert our expanded Dirac operator D into 
the Ginsparg- Wilson equation. To that purpose we de- 
fine: 



E = - D - 75-D75 
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and finding a solution D of the Ginsparg- Wilson equation 
corresponds to having E = 0. Computing the linear part 
of E is straightforward: When evaluating the sandwich 
75-D75 w e find that the terms of (Q) with an odd number 
of 7 P , i.e. vector- and pseudovector terms, pick up a mi- 
nus sign, while the other terms remain unchanged. Thus 
when adding the two linear terms we find that the terms 
with an odd number of 7^ cancel while the other terms 
pick up a factor of 2. 
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The next step is to compute the quadratic term 
75U75Z?. Here we have to multiply the various terms 
appearing in D. Each term is made of two parts, a 
generator of the Clifford algebra and a group of paths. 
The multiplication of two of these terms proceeds in two 
steps: Firstly the two elements of the Clifford algebra 
are multiplied giving again an element of the algebra. 
In the second step we have to multiply the paths of our 
two terms. This multiplication can be noted very conve- 
niently in our notation, where multiplication of two paths 
simply consists of writing the paths into one long path: 

< h, fa—In > x < hih-^n' >— < h>l2---ln:hih---ln' >- 
It is straightforward to establish this rule by translating 

back to the algebraic expression of our examples (0),(||) 

and performing the multiplication in this notation. It 

can happen that, after multiplying two paths, a hop in 

some direction li is immediately followed by its inverse 

—li. These two terms cancel each other and we find 

< h-.-h-i, Hi ~hi h+i---ln >=< h---H-i, h+i'-ln >• This 
rule is used to reduce all products of paths appearing 
in 75D75D to their true length. In a final step we de- 
compose the sums over the product terms into groups of 
paths related by the symmetries in the same way as we 
did above when constructing the most general ansatz for 
D. Adding the linear and quadratic terms of (^|) we end 
up with the following expansion for E: 
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Again the expansion of E is an infinite series and we dis- 
play only the leading groups of paths. The coefficients ef 



are now quadratic polynomials in the original coefficients 
Si,Vi,U, a,i and pi given by 

e\ = -2s x + sj + 8sl + 48sl + 8s 2 4 + 8v{ + 96v% + 8v% 

+A8tj + I92al + 38Apl ... , 
e 2 — — 2s 2 + 2si,s 2 + I2s 2 s 3 + 2s 2 s 4 + I2v x v 2 

+2v x v 3 ... , 

e s 3 = -2s 3 + 2s x s 3 + sj+ isj + 2s 3 s 4 + Av\ + 2v 2 v 3 ... , 
e| = -2s 4 + 2.S1S4 + s 2 + 6S3 - v\ - U\ - 2Aa\ 
-48p? ... , 

el = -S2V1 - As 3 v 2 - 2s 4 v 2 - s 3 w 3 - v 3 ti - 4v 2 h ... , 
e\ = -2t 1 + 2s 1 t 1 - 2s 4 ii - vf - Av\ - 2v 2 v 3 - U\ 
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e" = -s 2 ai + v 2 ti - v 3 pi ■■■ , 

e 2 = -v 2 h ... , 

e 3 = -s 2 ai - 2v 2 pi ... , 

e 4 = -2s 2 ai + 2w 2 <i - <iv 2 pi ... , 

e\ = —2p\ + 2s\pi — 2s 4 pi — 2v\a\ + t\ 



(7) 



Each coefficient ef is an infinite series of terms. 

The crucial step in our construction is to realize that 
the groups of paths together with the corresponding ele- 
ment of the Clifford algebra, appearing in the expansion 
(||), are linearly independent for arbitrary background 
gauge field. These groups of paths together with the 
products of 7 M have to be viewed as basis elements for the 
expansion of E. For a solution of the Ginsparg- Wilson 
equation we must have E = and thus all coefficients ef 
have to vanish simultaneously. Hence we have rewritten 
the problem of finding a solution of the Ginsparg- Wilson 
equation to the problem of solving the system of coupled 
quadratic equations (Q). 

Before we discuss solving (]?]) let us first analyze our 
Dirac operator D for trivial background gauge field. In 
this case we can compute the Fourier transform D(p) of 
D. For small momenta we have to implement the condi- 
tion D{p) — i p '+ 0(p 2 ) , which leads to two more equa- 
tions for the coefficients, the first one makes the constant 
term vanish and the second one requires the slope to be 
equal to 1. 



= si + 8s 2 + 48s 3 + 8s 4 .. 

1 = 2?7i + 24w 2 + 4w 3 ... . 



(8) 



These two equations have to be implemented as well, and 
together with (^) entirely describe Dirac operators which 
obey the Ginsparg- Wilson equation and have the correct 
behavior at small momenta. 

The quadratic system of equations (Q), (||) is com- 
pletely equivalent to the Ginsparg- Wilson equation. The 
known solutions of the Ginsparg- Wilson equation, i.e. the 
overlap operator as well as the perfect action, correspond 
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to two different solutions of (0) and (||). Many interest- 
ing questions, such as the existence of a connected con- 
tinuum of solutions of (R), (||) can be addressed in our 
new framework. 

In the remaining part of this letter we now explore the 
possibility of finding solutions of our system of quadratic 
equations numerically. It has to be stressed, that (0), 
(0) is a system of infinitely many quadratic equations for 
infinitely many coefficients. Thus in order to be able to 
apply numerical methods for solving these equations one 
has to introduce a cutoff, i.e. we set the coefficients for 
terms with longer paths to zero. This cutoff comes in 
a natural way for proper lattice Dirac operators, since 
in order to maintain universality of our lattice fermions, 
the coefficients Si,Vi,ti, a,i,pi have to decrease in size ex- 
ponentially as the length of the corresponding paths in- 
creases. 

In addition, the following qualitative argument (we use 
only the leading order, i.e. we ignore effects of back- 
folding paths) shows that the structure of the Ginsparg- 
Wilson equation strongly supports this exponential sup- 
pression of longer paths: When setting E = 0, the 
Ginsparg- Wilson equation (||) turns to D + 75-D75 = 
IbD^D. Short paths contained in the original D 
combine to longer paths when evaluating the product 
on the right-hand side. The linear term on the 
left-hand side requires also these longer paths with their 
coefficients to be contained in D + 75-D75. The coeffi- 
cients of the combined longer paths are products of the 
coefficients of the short paths, and since the coefficients 
of the shortest paths are already smaller than 1 we find 
an exponential suppression of the coefficients for longer 
paths. This behavior was confirmed numerically |]||. 

Thus the physical requirement of universality as well 
as the structure of the Ginsparg- Wilson equation imply 
the coefficients to decrease exponentially in size as the 
length of the paths in D increases. We now consider a 
truncated expansion of D which contains only terms with 
paths up to a certain length. Introducing this cutoff is 
equivalent to what is done when computing numerically 
the fixed point action using a finite parametrization for 
the Dirac operator. The resulting Dirac operator will 
only be an approximate solution of the Ginsparg- Wilson 
equation, but the error, i.e. the norm of E, can be made 
arbitrarily small by including more terms in D. 

After one has introduced the cutoff discussed above, 
e.g. one includes only terms with a maximal length of 3, 
the system (Q) , (||) turns into a finite system. Each coeffi- 
cient ef in the system (fjj) as well as the two equations (||) 
now consist of only finitely many terms. We numerically 
solve the equations corresponding to the leading terms 
in the expansion of E, and so determine the coefficients 
Sj, Vi, ti, a>i,pi in the expansion (^) of D. 

This approach has been implemented and carefully 
tested in two dimensions || . We found that adding only 
a few simple terms to the Dirac operator improves the 



chiral properties tremendously, at only a small increase 
of the cost of a numerical treatment. The resulting Dirac 
operator was tested in a dynamical simulation of the 2- 
flavor Schwinger model and it was found that the mass of 
the 7r-particle in this model can be brought down, close 
to the values obtained with the perfect action and the 
overlap operator. 

We are currently implementing the quadratic system in 
4 dimensions and have begun to construct approximate 
solutions of the Ginsparg- Wilson equation along the lines 
outlined in this letter. Preliminary results are very en- 
couraging, in particular we find a considerable reduction 
of the additive fermion mass renormalization. These re- 
sults will be presented elsewhere jj). 
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